Order-to- disorder transition in ring-shaped colloidal stains 
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A colloidal dispersion droplet evaporating from a surface, such as a drying coffee drop, leaves a 
distinct ring-shaped stain. Although this mechanism is frequently used for particle self-assembly, the 
conditions for crystallization have remained unclear. Our experiments with monodisperse colloidal 
particles reveal a structural transition in the stain, from ordered crystals to disordered packings. 
We show that this sharp transition originates from a temporal singularity of the flow velocity inside 
the evaporating droplet at the end of its life. When the deposition speed is low, particles have time 
to arrange by Brownian motion, while at the end, high-speed particles are jammed into a disordered 
phase. 



Ordered arrays of colloidal particles present important 
characteristics for fields as photonics [ ] and biotechnol- 
ogy [ ] . A relatively simple approach to deposit particles 
onto a substrate is by evaporation of a colloidal dispersion 
droplet. In droplets with pinned contact lines, evapora- 
tion gives rise to the so-called coffee-stain effect [^-s]: a 
ring-like stain remains on the substrate once the liquid 
has evaporated. The mechanism behind this ring-stain 
formation is explained in the pioneering work of Dee- 
gan et al. [ ]: In an evaporating drop with an immobile 
contact line, a capillary flow is generated to replenish 
the liquid that has evaporated from the edges. This flow 
drags particles towards the contact line, forming the ring- 
shaped stain in particle suspensions such as coffee. 

This deposition mechanism has been successfully used 
to generate colloidal crystals in three dimensions [ ] . Fur- 
thermore, techniques have been developed to control the 
patterns by permitting evaporation through specially de- 
signed masks [ ] or by controlling the wettability [ , ] . 
Also large-scale colloidal crystals have been observed for 
particle sizes below 100 nm [ ], but some displayed in- 
stabilities that gave rise to concentric rings [ ], stripes, 
chevrons [ ], and cracks [ ]. At present, there is no way 
of predicting beforehand the structure of the deposit re- 
maining after evaporation, due to a lack of understanding 
of the basic mechanisms. However, the processes driv- 
ing colloidal crystallization are crucial for controlling the 
mechanical and the optical properties of the synthesized 
materials. 

In this Letter, we reveal the conditions leading to 
colloidal crystallization in a paradigmatic system: a 3- 
/iL sessile water droplet containing colloidal particles 
evaporating from a smooth glass slide (Fig. la). The 
colloids consist of almost monodisperse red-fluorescent 
polystyrene particles with particle diameters in the range 
of 0.5 to 2 jivd. The droplet is gently deposited onto a 
glass, and is left to evaporate under controlled temper- 
ature and humidity. The contact line remains pinned 
during the entire experiment. Our aim is to analyze 
structure of the particle arrangement in the remaining 




FIG. 1. Order-to-disorder transition in the particle stain left 
by an evaporating drop, (a) A 3-/iL sessile water droplet evap- 
orates from a glass substrate, (b) Ring-shaped stain of red 
particles (the coffee stain) left on the substrate after evapora- 
tion, (c) A close-up by an optical microscope of the bottom 
layer of the stain, taken from the white square in (b), shows 
that the outermost lines of the stain (left) have an ordered, 
crystalline structure. Towards the center of the drop (right), 
a transition to a disordered particle arrangement is observed. 
(d) A top view of the ring stain, taken from the red square 
in (c) with the Scanning Electron Microscope (SEM), shows 
that the first lines of particles (left) are arranged in a hexag- 
onal array, while the following lines (brighter in the image) 
are arranged as square, followed by again hexagonal array. 



coffee stain. It turns out that the dynamics of particle 
deposition is crucial for the understanding of the final 
structure. Therefore, the entire evaporative process is 
filmed simultaneously from the side and from below [ ] . 
This allows to follow the evolution of the drop geome- 
try, the formation of the deposit, and the dynamics of 
particle migration. The particles are filmed from below 
with a CMOS camera (PCO imaging) connected to an in- 
verted microscope. An integrated standard halogen light 
source is used to illuminate the particles; the emitted 



light is filtered before it reaches the camera. The camera 
is operated at a typical frame rate of 1 frame per second, 
which is sufficient to capture the movement of the parti- 
cles. To determine the velocity field in the measurement 
area from the captured images, a custom-made MAT- 
LAB micro-Particle Image Velocimetry (/iPIV) code was 
used. The side-view images are used to calculate contact 
angle, volume, radius, and height of the drop at every 
instant with a custom-made MATLAB code. Side-view 
images are taken with a long distance microscope and a 
CCD camera (Lumenera Corp.), synchronized with the 
visualization of the particles by a common trigger. 

When the droplet has completely evaporated, most of 
the colloidal particles have aggregated into a ring-shaped 
coffee stain (Fig. lb). Zooming in on the stain, we see 
that the particle arrangement is not homogeneous (Fig. 
Ic): there is a remarkable transition from a crystalline 
arrangement to a disordered phase. A top view of the 
stain (Fig. Id) shows a sequence of particle arrange- 
ments within the ordered phase, starting from hexagonal 
packing, followed by square packing, and again hexag- 
onal packing. The particle arrangement is analyzed in 
more detail in Fig. 2. Close to the contact line, we 
find a square arrangement of particles (Fig. 2bl), fol- 
lowed by a hexagonal array (Fig. 2b2), and finally a 
disordered phase (Fig. 2b3). The same trend has been 
observed for different particle sizes and concentrations. 
A well-known method to quantitatively assess the differ- 
ent packing structures, is to calculate the Voronoi areas 
[17] of the particles, as shown in Fig. 2b. In Fig. 2c, we 
identify the Voronoi areas corresponding to the different 
cell arrangements: the disordered phase exhibits a larger 
average cell area as well as a wider dispersion, compared 
to the ordered phase. This reveals the sharp transition 
from order to disorder. 

The order-to-disorder transition can be traced back to 
the hydrodynamics within the drop. Fig. 3a shows the 
radial velocity measured by /iPIV of the layer of particles 
closest to the glass slide, in a region close to the contact 
line. The particle velocity increases dramatically in the 
last moments of the droplet 's life. We refer to this sud- 
den change in speed as "rush hour". The particles that 
arrive early, at a low deposition speed, form an ordered 
(square or hexagonal) structure. By contrast, particles 
that arrive during rush hour have a high speed and form 
a jammed, disordered phase. 

The rush-hour behavior can be explained by simple 
mass-balance considerations. The volume flow towards 
the contact line inside the drop is driven by the evapo- 
ration from the drop surface ['^]. It turns out that the 
rate of evaporation is approximately constant over time 
[6]. To replenish this evaporated liquid, a continuous 
volume flow towards the contact line is generated inside 
the drop. However, the drop height is vanishing during 
evaporation. This can be characterized by a contact an- 
gle 0{t) -^ at the final stages of evaporation (see Fig. 
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FIG. 2. Analysis of the particle ordering in the stain, (a) 
Overview of the different patterns observed in the ring stain. 
The position of the contact line is on the left side of the figure 
(the first layers of particles are not shown), (bl-3) A close-up 
reveals the sequence of patterns in the stain: square packing 
close to the contact line (bl), followed by hexagonal packing 
(b2), and finally disordered packing (b3). The Voronoi cells 
belonging to the particles are shown in blue, (c) The area 
of the Voronoi cells plotted versus the distance x from the 
contact line. The order-to-disorder transition at x = Xc is 
determined by the increase in the mean Voronoi area and 
the dispersion, due to the lack of order. This transition is 
indicated by the red dotted line. 



3). Hence, the same amount of liquid has to be squeezed 
through an area which is vanishing, inducing a diverging 
radial velocity. Close to the contact line, the dependence 
of the height- averaged radial velocity u on time t and 
distance from the drop center r, is given by [ ., v.] 



u{r,t) = 



D* 



0{t) ^R{R-ry 



(1) 



see l'^^] for details. Here, D* = 2^/2Dya^c/{7rp) repre- 
sents the driving of the flow, with Dya = 24 x 10~^ m^/s 



3 
"w" 2 

£ 

3 1 

0! 

-1 



1 l/^t^ ^-^^^ • • • •• • • • ^ 


\ 0(t 


i 


n 


^ 


J 




p 


..0^-^ 



500 1000 I, 2000 

evaporation time(s) 

FIG. 3. Origin of the order- to- disorder transition. Plot of the 
radial particle velocity versus time, measured at a distance of 
174 /im from the contact line. The fluid velocity predicted 
by the model (1) is plotted as a solid line. Both model and 
experiment show the dramatic velocity increase at the end 
of the droplet 's life (the rush hour). The inset shows the 
direction of the radial velocity in the droplet. The time tc at 
which the order-to-disorder transition occurs is determined 
from the videos; the corresponding critical velocity Uc is read 
from this graph (indicated by the black line in the figure). 
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FIG. 4. The experimentally determined critical velocity Uc 
versus its theoretical prediction (3). 



the diffusion constant for vapor in air, Ac = 1.2 x 10~^ 
kg/m^ the vapor concentration difference between the 
drop surface and the surroundings, p = 998 kg/m^ 
the liquid density, and R the drop radius. The spa- 

-I /o 

tial {R — r) ' -divergence originates from the diverging 
evaporative flux from the drop surface [3] . Although this 
spatial divergence generated a lot of interest ['^-^], it is 
clear from (1) that there is also a temporal divergence be- 
cause 0{t) vanishes. This temporal divergence provides 
the key to the rush-hour behavior: the radial velocity 



blows up towards the end of the droplet 's life. 

The experimental velocity measurements are per- 
formed at a distance of 2 fim above the glass slide, and 
therefore differ from the theoretical height- averaged ve- 
locity. To quantitatively compare the model with the 
experimental data, a description of the velocity profile 
is required. Given the low velocities measured inside 
the evaporating drop (low Reynolds number) the inertial 
terms in the Navier-Stokes equation can be neglected. 
Since the droplet height is much smaller than its radius, 
the thin-film (lubrication) approximation [^ ^] can be used 
to describe the velocity field. Hence, the velocity u{r^ z, t) 
[16] is given by 



u{r,z,t) 



^^n{r,t)(h{r,t)z-\z-], (2) 



with h{r,t) = 0{t){R'^ -r^) /{2R), the local droplet 
height. As shown in Fig. 3, this model quantitatively 
describes the rush-hour behavior. For the plot of (2) 
we took R — r = 174 /am. (the distance from the con- 
tact line where the /i-PIV measurements are performed), 
z = 2 /im and t^ the droplet 's life time determined from 
the experiment. This means there are no adjustable pa- 
rameters in the comparison between the model and the 
experiments. 

We show that the structural change in the stain can be 
entirely attributed to the rush-hour behavior. The par- 
ticles that form the ordered phase arrive at the contact 
line at an early time, when the deposition rate is still low. 
On the other hand, the particles that belong to the disor- 
dered phase arrive during rush hour, at a high deposition 
rate. The slow, early arriving particles have time to ar- 
range themselves within the stain by Brownian motion, 
while those coming in rush hour have no time to find 
an appropriate spot and are quenched into a disordered 
phase. We can quantify this transition by comparing the 
time scale on which particles arrange by diffusion to the 
time scale of particle deposition. The diffusive time scale 
is to = dp/D^ where D is the diffusivity of the particles 
in liquid, calculated by the Einstein relation. The hydro- 
dynamic time scale is th = L/u (where L is the typical 
distance between the particles, dependent on their con- 



centration in the solution np{t) as L 



Til 



-1/3 



, and u is 



the particle velocity). When t^ and t^ become of the 
same order, a transition is expected. We therefore define 
the critical transition velocity as 



LD/dl. 



(3) 



We validated this criterion experimentally for various 
particle concentrations and sizes. First, we determined 
the time tc when the ring-shaped stain undergoes the 
transition (the position where this occurs is Xc, as in- 
dicated by the red line in Fig. 2c). Next, we used the 
velocity-time plots (Fig. 3) to find the experimental value 



ofu{tc) = Uc. As demonstrated in Fig. 4, the experimen- 
tal critical velocity perfectly follows the proposed scaling. 
The agreement suggests that the pref actor in (3) is actu- 
ally of order unity. Our model also explains many results 
described in the literature [^ '^-^ ^], where large-scale crys- 
talline deposits were observed for particle sizes ranging 
from 6 to almost 50 nm, in hindsight. In this range of 
particle sizes, the critical velocity required for a disor- 
dered phase (3) is unreachable within the droplet, which 
explains the absence of a disordered phase. 

Apart from the order-to-disorder transition, there are 
also transitions within the ordered phase, from square to 
hexagonal packing and vice versa (Fig. Id). In an un- 
confined system, the most efficient packing lattice is the 
Hexagonal Close Packing (HCP), as conjectured by Ke- 
pler [ ]. However, in the evaporating drop the particles 
are confined in a wedge formed by the glass slide on one 
side, and the liquid-air interface with contact angle 6 on 
the other side. It was shown by Pieranski et al. [20, 2i] 
that such a confinement indeed leads to a sequence of 
hexagonal and square packing structures, dependent on 
the most efficient packing for the available space. In- 
deed, when a new layer is formed, the confinement by 
the wedge favors the square packing; see Fig 1. Away 
from the step edge, the available space increases and al- 
lows for the denser, hexagonal packing structure. We 
observe a sequence of this type of transitions, each time 
a new layer is added. Similar transitions were also ob- 
served by Abkarian et al. [22] in a wedge formed by a 
moving contact line. 

In conclusion, in this Letter we revealed the mechanism 
triggering crystallization in an evaporating colloidal so- 
lution drop with a pinned contact line. The temporal 
divergence (rush hour) observed in the velocity field in- 
side the drop is responsible for the lack of order in the 
ring-shaped stains: when the deposition rate is low, par- 
ticles have time to arrange by Brownian motion and form 
an ordered phase. During the rush hour, there is no time 
for such re-arrangement and the particles form a disor- 
dered phase. Such effect resembles the popular computer 
game Tetris, where at the beginning of the game slowly 
falling objects are easily arranged in regular structures, 
while as the game evolves, the fast falling objects become 
jammed and disordered. Therefore, the key to generate 
large-scale colloidal crystals is to avoid the rush-hour be- 
havior, i.e. to avoid small contact angles during evapora- 
tion. This can be achieved for example by modifying the 
hydrophobicity of the substrate or by continuously re- 
plenishing the evaporated liquid. Large-scale crystalline 
arrays can also be obtained by tuning the particle size, 
since nanometer-sized particles diffuse much more effi- 
ciently. As a rule of thumb, freely-evaporating colloidal 



solutions with particles within the nanometric range will 
naturally crystallize, micrometer-sized particles will show 
both ordered and disordered phases, while larger parti- 
cles lead to disordered stains 
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